ESD  TR- 67-463 
ESTO  FILE  COPlf 


11—6 


ESD-TR  67-463 


esd  accession  list 

tSTI  Call  No._  AL  5  789? 

Copy  NO. 

TRAVEL -TIME  ANOMALIES 
AT  LASA 


ESD  RECORD  COPY 

RETURN  TO 

SCIENTIFIC  &  TECHNICAL  INFORMATION  DIVISION 
(£STty,  Bl'ILDWG  1VL1 


7  JULY  1967 


REPORT  No.  L L  —  6 


Prepared  for 


LINCOLN  LABORATORIES 

MASSACHUSETTS  INSTITUTE  of  TECHNOLOGY 


EARTH  SCIENCES 

A  TEI.EDYNE  COMPANY 


AD  653  Hn. 


Hub  document  lias  been  approved  for  public  release  "and  "Baioj 
its  distribution  is  unlimited. 


59 


TRAVEL -TIME  ANOMALIES 
AT  LASA 


7  JULY  1967 


REPORT  No.  L  L  —  6 


Prepared  for 


L  I  NCOLN  L  A  BOR  A  TOR  I  ES 
MASSACHUSETTS  INSTITUTE  of  TECHNOLOGY 

Purchase  Order  No*  BB-246 

PRIME  CONTRACT  AF  19(628)-5167 


EARTH  SCIENCES 


ATFI  EOYNE  (  OMPANY 


TABLE  OF  CONTENTS 


Page  No 


FOREWORD 

ABSTRACT 

1.  INTRODUCTION  1 

2  .  PROCEDURE  1 

2.1  Data  Processing  1 

2.2  Error  Detection  3 

2.3  Computer  Programs  5 

3.  CONCLUSIONS  AND  FUTURE  PLANS  5 

ILLUSTRATIONS 
APPENDIX  A 
APPENDIX  B 


FOREWORD 


The  work  documented  in  this  report  is  part  of  a  study  of  the  travel¬ 
time  and  amplitude  anomalies  observed  at  the  Large  Aperture  Seismic  Array 
(LASA)  in  Montana.  This  report  supersedes  Report  No.  LL-1 ,  "Travel-  Time 
Anomalies  at  LASA"  in  that  additional  data  and  data  quality  control  have 
been  used  to  obtain  the  results  presented  in  this  report. 

This  report  was  written  by  D.  E.  Frankowski  and  F.  A.  Klappenberger . 
Assistance  in  the  collection  of  data  was  provided  by  A.  L.  Kurtz,  V.  R. 
McLamore ,  R.  D.  Mierley,  P.  A.  Santiago,  H.  N.  Johnson,  F.  J.  Whited, 
and  D.  A.  Beschen.  Dr.  E.  F.  Chiburis  served  as  a  consultant.  The  pro¬ 
ject  director  was  Dr.  P.  W.  Broome. 


ABSTRACT 


World-wide  average  travel-time  curves  are  inadequate  to  permit 
optimum  phasing  of  teleseismic  signals  Phase  alignment  corrections 
in  the  form  of  station  relative  travel-time  anomalies  of  LASA,  based 
on  approximately  570  events,  are  presented . 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 


1. 


INTRODUCTION 


World-wide  average  travel-time  curves  are  known  to  be  inadequate 
to  permit  optimum  phasing  of  teleseismic  signals  This  is  not  a  ser¬ 
ious  problem  for  events  of  large  recorded  earth  motion  since  it  is 
possible  to  adjust  a  preliminary  alignment  using  measurements  taken 
from  the  signal  itself.  For  events  where  the  signal  is  obscured  by 
noise,  however,  these  corrections  to  the  preliminary  alignment  must  be 
known  a  priori.  This  report  presents  time-shift  corrections  in  the 
form  of  station  relative  travel-time  anomalies  for  the  Montana  LASA. 
These  corrections  are  based  upon  travel-time  observations  at  the  center 
seismometers  of  all  21  subarrays  for  approximately  570  events  (see 
Figure  1) . 


2 .  PROCEDURE 

2 . 1  Data  Processing 

Figure  2a  illustrates  the  manner  in  which  the  initial  P-wave  from 
an  event  travels  from  the  source  to  LASA.  Ideally,  the  wave  front 
passes  across  LASA,  as  shown  in  Figure  2b.  Because  of  travel-time 
anomalies,  however,  this  idealized  situation  is  not  observed. 

Figure  2c  shows  a  typical  record  of  an  event  as  seen  on  the 
16mm  Develocorder  film.  Two  such  films  record  the  earth  motion  for 
all  21  center  seismometers.  Necessary  corrections  to  the  JB  travel¬ 
time  tables  are  established  by  comparing  these  arrival  times  with  the 

predicted  times.  This  establishes  corrections  for  a  particular  event 
1 

region . 

Data  recorded  for  one  event  is  shown  in  Figure  3a.  This  data  is 
keypunched  as  indicated  in  Figure  3b.  Then,  using  the  definition  of 
the  travel-time  anomaly 


Ot  =  (T  - 
I  I 


pr 


-  (To  - 


V 


1.  Chiburis  and  Dean,  "Teleseismic  Signal  Alignment  at  the  Tonto 
Forest  Extended  Array,"  SDL  Report  No.  125,  15  October  1965. 
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where  T  and  P  are  observed  and  predicted  travel-times  at  station  I 
and  reference  station  0,  respectively,  the  anomalies  are  computed 
using  AO-10  as  a  reference  as  shown  in  Figure  3c. 

In  all  cases,  the  USC&GS  source  data  and  JB  travel-time  curves 
are  used  as  a  standard.  Any  bias  errors  in  USC&GS  data  will  thus  be 
reflected  as  correction  errors.  The  P-phase  arrival  times  at  each 
center  are  read  from  the  16mm  film  at  a  magnification  of  20  on  a 
projected  image  view  screen.  The  timing  point  of  each  arrival  is 
defined  as:  "the  earliest  motion  within  the  first  quarter  cycle  which 
can  be  read  most  accurately  and  consistently  on  all  channels." 

Since  some  channels  may  have  poor  traces  the  readings  are  all 
qualified  by  a  quality  grade  as  follows. 

Quality  Grade  1  =  a  perfect  trace 

Grade  2  =  a  trace  that  is  easily  read 

Grade  3  =  a  trace  on  which  the  signal  is  seen  but 
is  not  distinct 

Grade  4  =  a  trace  on  which  the  signal  cannot  be  seen. 

The  object  of  this  process  is  to  obtain  corrections  which  can  be 
defined  as  a  mathematically  simple  function  of  distance  and  azimuth 
from  LASA .  To  achieve  this,  events  were  grouped  into  cells  by  distance 
and  azimuth  relative  to  LASA  and  the  average  relative  travel-time  anom¬ 
aly  of  all  events  in  a  cell  were  then  calculated  for  each  subarray. 

The  travel-time  anomaly  was  assumed  constant  within  each  cell  at  this 
point  in  the  analysis.  The  standard  deviation  and  number  of  events 
included  in  the  cell  population  were  also  calculated  (a  sample  of  this 
result  is  shown  in  Figure  4).  The  complete  anomaly  data  is  shown  in 
Figures  6  through  25. 

One  would  expect  large  standard  deviations  in  a  cell  population  if 
either  the  anomaly  were  not  nearly  constant  within  a  cell,  or  if  ab¬ 
normal  errors  existed  in  the  data  (The  results  in  this  report  were 
obtained  by  using  all  except  Grade  4  data) . 
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2.2  Error  Detection 


The  desired  result  has  been  to  obtain  average  travel-time  anomalies 
to  within  0.1  second  or  better.  Using  the  data  from  Figures  6  to  25, 
a  statistical  study  was  undertaken  so  that  outliers  could  be  defined 
and  re-examined.  The  error  model  used  in  this  study  was  based  on  the 
classical  notion  of  experimental  error  due  to  inherent  observation 
variations  The  error  model  proposed  in  this  case  is  the  result  of  two 
main  causes.  The  first  component  of  the  model  is  due  to  the  error  in 
reading  the  film  traces.  The  second  is  due  to  epicenter  location  errors. 
We  assume  these  errors  to  be  normally  distributed*'  and  independent  of 
one  another.  Thus,  the  total  reading  error  Z  is  assumed  to  be  such 
that 
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where  U  and  U  ,  c  and  a  are  the  means  and  variances  of  the  reading 
R  L  R  L 

and  location  errors  respectively. 

The  mean  due  to  reading  is  assumed  to  be  zero.  Then,  since  all 
readings  are  referred  to  the  USC&GS  locations,  we  further  assume  that 


U  is  zero.  The  effect  of  this  is  to  correct  LASA  so  that  it  agrees 
with  USC&GS • 


An  estimate  S  ,  of  °  ,  the  variance  of  the  entire  population,  is 
computed  in  each  cell.  By  pooling  these  estimates  one  obtains  an  im¬ 
proved  estimate  which  is  very  good  due  to  the  large  number  of  cells. 
The  procedure  is  as  follows:  The  estimates  in  each  cell  are  entered 
on  a  histogram,  Figure  5,  as  a  function  of  the  number  of  events  in 
that  cell . 

A  pooled  estimate  of  the  population  variance  is  then 
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*The  statistical  methods  used  are  not  especially  sensitive  to  the 
accuracy  of  this  assumption. 


where  f  is  the  number  of  estimates  in  the 
NO 

by  j  and  with  N  events  in  each  of  the  cells. 


class  interval  indexed 

2 

The  value  of  S  is 


D 

equal  to  the  center  value  of  the  variance  in  the  class  of  interval  j . 

2 

The  estimate  for  a  is  sufficiently  good  so  that  we  may  compute 

2 

100(l-a)%  confidence  intervals  for  S  and  S  „ 


2  -2 

We  set  a  =  S  , 

2  ^ 

and  o  =  (Sz)  . 

The  confidence  intervals  are  then  defined  from  the  relation. 
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where  f  =  number  of  degrees  of  freedom  and  Xz 


(a/2) ,f 


and  X  , ,  „  .  _  are  found  in  tables 

(l-a/2) , f 


After  these  confidence  limits  are  established,  we  then  proceeded 
to  question  all  outliers.  Since  °  is  estimated  from  the  population  and 
it  is  the  population  which  are  questioning,  this  is  a  bootstrap  oper¬ 
ation.  If,  for  N  events,  a  was  computed  which  fell  within  the  95% 

2  2 
limits,  the  S  was  accepted  as  a  reasonable  estimate  for  °  .  If,  how- 

2 

ever,  the  computed  S  fell  outside  the  confidence  intervals,  the  esti¬ 
mate  was  questioned  and  the  time  anomalies  involved  in  computing  the 
estimate  were  re-examined  for  errors  beyond  the  expected  film  reading 
and  epicenter  location  errors.  In  this  way,  about  1%  of  the  readings 
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were  found  to  be  in  error.  These  errors  have  been  removed  from  the 
data  in  Figures  5  through  26. 

2 . 3  Comput er  Programs 

Two  programs  were  used  in  this  study.  One  computed  the  anomalies 
using  the  film  readings  and  the  USC&GS  event  data.  The  second  computed 
average  anomalies  in  specified  range  and  azimuth  cells  given  the  anomaly 
data.  The  second,  in  addition,  computed  estimates  of  the  error  variance 
in  each  cell.  Program  write-ups  are  included  as  Appendices  A  and  B. 

3 .  CONCLUSIONS  ANT)  FUTURE  PLANS 

Travel-time  anomaly  data  for  geographic  regions  of  high  seismicity 
can  be  taken  from  Figures  6  through  25  or  can  be  supplied  in  IBM  card 
form  by  region  or  on  an  event-by-event  basis.  This  data  may  contain 
some  errors.  We  are  continuing  to  check  the  data  for  obvious  blunders. 
Results  at  all  times  where  there  are  few  events  must  be  considered  as 
tentative  since  the  only  check  we  have  on  USC&GS  errors  is  use  of  re¬ 
peated  events. 

As  more  events  are  added,  it  becomes  easier  to  eliminate  both 
reading  and  basic  data  errors.  We  will  continue  to  investigate  those 
cases  where  the  data  in  a  cell  exhibits  high  variability  and  where  no 
obvious  error  has  been  found.  Since  the  mean  travel-time  anomaly 
varies  from  cell  to  cell  it  is  not  unreasonable  to  expect  that  one 
subdivision  of  a  cell  may  exhibit  a  mean  anomaly  clearly  dissimilar 
to  that  of  another  portion  of  that  same  cell. 


5 
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Figure  1.  Event  location  chart 
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Figure  2c.  Typical  recording  showing  event  arrival 
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Figure  3c.  Computer  output  listing  computed  Figure  3b.  Observed  data 

anomalies  for  one  event  in  keypunched  form 
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Figure  4.  Sample  output 


34 


C 

o 

g 

> 

<D 

TD 


X) 

%— 

O 

TJ 

C 

o 


<n 


32 


30 


20 


26 


24 


.22 


.20 


.18 


.16 


.14 


.12 


.10 


08 


.06 


04 


.02 


.00 


LEGEND 

SIZE  OF  BLOCK  SYMBOLIC 
OF  NUMBER  OF  POINTS 
AT  THE  INTERSECTION 


|  «  SIX  DATA  POINTS 


✓ 


2  4  6  8  10  12  14  16  18  20  22  24  26  28  30 

N  (number  of  events) 


Figure  5 


Distribution  of  standard  deviation 
computed  from  N  events 


3  000  4000  5000  6000  7000  8000  9000 


o 

ro 


8 


I 


o  o 

id  <r> 

- - *-  1 


o 

CM 


O 

if) 


o 

00 


o 

CM 


o 

O 

O 

O 

O 

N 

O 

ro 

CD 

CM 

CM 

ro 

ro 

ro 

■  O  T  O 


"T" 

O 


y- 


\ 


x 


X? 


O 

CD 


O 

(T> 


O 

CVJ 


-?4- 


X 


o 

.8 


\ 


®  PQ 


<D 

E 

o 


O 

o  § 
o  < 
o 

O 

cr 


LU 

o 

o  < 

O  I- 

2  </> 

Q 


rd 

S-l 

S4 

rd 

rQ 

D 

c/i 


vD 

QJ 

Jh 

3 

On 


O 

■8 

ro 


o 

ID 


o 

GO 


o 

CM 


O 

P 

o 

O 

<3* 

r* 

o 

ro 

CM 

CVJ 

ro 

ro 

O 

8 


(saej63p)  QI0V  01  SAIlVISd  HIOWIZV 


o 

ro 


O 

<D 


I 


I 


o 

o_ 

o 

0> 


o 

o_ 

o 

u> 


o 

J- 


o 

CJ 


o 

m 


o 

oo 


o 

CJ 


o 

o 

o 

o 

O 

h- 

o 

ro 

CD 

CJ 

OJ 

ro 

ro 

ro 

o 

ro 


“P" 

O 

<D 


V- 


\ 


X 


X; 


y 


X 


o 

■8 


\ 


to 

k_ 

<t> 

CN 

(D 

CQ 

E 

o 

>i 

rtf 

w 

o 

rtf 

< 

rQ 

3 

Cfl 

o 

cr 

0 

u_ 

LU 

<U 

o 

U 

3 

< 

Cn 

H 

rH 

C/) 

Cu 

Q 

O 

in 


O 

00 


o 

OJ 


o 

p 

o 

O 

F- 

o 

ro 

CJ 

CM 

ro 

ro 

o 

s 


(saajBap)  010V  01  3AllV~)3a  H1DIMIZV 


3  000 _ 4000 _ 5000  _ 8000  7000 _  8000 _ 9000 _ 10,900 _ 11,900 


o 

ro 


O 

<0 

JL 


O 


o 

c\j 


o 


o 

00 


o 

PJ 


o 

o 

o 

o 

O 

r- 

o 

ro 

iO 

PJ 

pj 

ro 

ro 

K) 

4 


¥- 


\ 


X 


X 


7 


ss: 


.y 


X 


o 

1-8 

o' 


o 

■8 

CD 


\ 


<D 

ro 

00 

CD 

CQ 

E 

o 

>i 

— 

rd 

DC 

Jh 

O 

rd 

<3 

rQ 

O 

o 

< 

3 

O 

cn 

5 

O 

cr 

° 

Li_ 

00 

LU 

0) 

O 

V-i 

z 

a 

O 

< 

O 

O 

(D 

H 

00 

rH 

Pu 

o 

o 

O 

O 

O 

o 

o 

o 

o 

o 

o 

O 

ro 

CD 

CD 

PJ 

ID 

go 

o 

to 

CVJ 

pj 

pj 

ro 

to 

o 

8 


(S99J69P)  010V  01  3 AllV13d  HIOWIZV 


3  000  4000  5000  6000  7000  8000  9000  10,000  11,000 


o 

ro 


O 

iD 


o  6  3  0 


O 

(7) 

-L 


O 

Cd 


o 

m 


o 

co 


o 

Cd 

_L_ 


o 

o 

O 

O 

O 

<*• 

h* 

O 

ro 

(£> 

CM 

Cd 

ro 

ro 

ro 

-A 


o 

ro 


o 

to 


V- 


X 


V 


X 


-74- 


X 


o 

1-8 

o 


o 

■8 


\ 


to 

k_ 

<D 

<D 

E 

o 


O 

5> 


O 

CT 


CQ 

>i 

U 

u 

rtf 

3 

cn 


Ll 


LU 

O 

i< 

S  <2 

Q 


Q) 

M 

3 

tn 

rH 

PU 


o 

in 


o 

00 


o 

Cd 


O 

p 

O 

o 

*f 

r~ 

O 

ro 

Cd 

Cd 

ro 

ro 

O 

s 


(S99J63P)  010V  01  3AI1V33U  HlflWIZV 


3  000  4000 _ 5000  6000 _ 7000 _ 8000 _ 9000 _ 10,900 _ H,000 


o 

O  ro 

L_— - 


o 

t£> 


i 

o 

-J 


O 

CT> 

mJLm 


o 

CVJ 


O 

«■) 


O 

00 


o 

CVJ 


o 

o 

o 

o 

O 

h- 

o 

ro 

u> 

CVJ 

CVJ 

ro 

ro 

ro 

-rt-3  M 


-  X 


7| 


X 


■V 


X 


\ 


iA 

k_ 

<D 

CD 

£ 

o 


g 

o  ® 
o  < 
o 

o 

cn 


LU 

o 

§5 

S  « 

Q 


O 

O 

o 

o 

o 

o 

o 

o 

o 

R 

O 

O 

to 

<0 

0) 

(VI 

in 

00 

CVJ 

<3- 

O 

ro 

CVJ 

CVJ 

ro 

ro 

O 

s 


(saajBap)  010V  01  3AllV13d  H10IMIZV 


U 

>1 

03 

U 

U 

03 

-Q 

2 

co 


CD 

w 

3 

tn 

r-} 

Pn 


o 

ro 


O 

CO 


o 

o> 


o 

OJ 


o 

m 


o 

05 


O 

CM 


o 

o 

O 

O 

O 

'f 

h- 

o 

ro 

CO 

OJ 

0J 

ro 

ro 

ro 

8. 

o.- 


O 

o_ 

o 

to 


-A 


V- 


\ 


v 


v 


'V 


o 
■  8 


O 

■8 


\ 


O 

VS 

ro 


O 

o 

O 

o 

o 

O 

O 

6 

O 

R 

CM 

O 

O 

ro 

CO 

OJ 

in 

GO 

CM 

OJ 

O 

ro 

to 

ro 

o 

s 


(saejBap)  OI0V  01  3AllV13y  H10WIZV 


to 

k_ 

CD 

Cm 

CD 

O 

£ 

o 

>i 

03 

Vj 

w> 

O 

03 

<s 

< 

£ 

P 

U) 

o 

o: 

u_ 

rH 

*H 

UJ 

o 

0) 

~2L 

H 

< 

D 

H 

tn 

CO 

H 

O 

Cm 

000*11  0o6‘0l  0006  0008  oool  0009  OOOS OOOb  oooT 


o 

jo 


o 

<0 


O 

cn 

_L_ 


O 

OJ 


o 

m 


o 

CD 


O 

OJ 


o 

o 

o 

o 

o 

h- 

o 

JO 

<0 

0J 

0J 

JO 

JO 

JO 

4 


o 

JO 


~~r~ 

o 

<0 


¥- 


\ 


V 


7^ 


o 

in 


o 

CD 


O 

0J 


T" 

O 

<*• 


"T" 

R 

0J 


X 


o 

-8 
a ) 


\ 


o 

o 

o 

JO 

JO 

JO 

co 

w 

<D 

a> 

E 

o 


O 

<S> 


O 

cr 


UJ 

a 


co 


00 

u 

>i 

u 

u 

03 

D 

CO 


(N 

rH 

a> 

D 

Cn 


Q 


O 

_o 

8 


o 

YS 

JO 


O 

s 


(saajBap)  010V  01  3AllVH3d  HIOWIZV 


3  0,00  4000  5000  6000  7000  8000  9000  10,000  11,000 

■  ■  1  f  f  I  -  i  |T  i  ■  .  i  i  I—  .  ■  ■  T  ,  ,i 


o 

ro 


O 

to 


O 

CD 

-L. 


O 

CM 


O 


o 

co 


o 

CM 


o 

O 

O 

O 

O 

h- 

o 

ro 

to 

CM 

CM 

rO 

ro 

ro 

y- 


\ 


■v 


V 


ss: 


V: 


X 


\ 


o 

■8 

CD 


if) 

k_ 

<0 

< 0 

u 

E 

o 

>1 

ffl 

U 

' — 

u 

O 

rt 

S> 

< 

rQ 

D 

CO 

o 

cr 

# 

u. 

ro 

t-H 

UJ 

o 

<U 

z 

U 

< 

P 

1- 

tn 

c n 

•  H 

o 

tn 

.8 

8 


o 

o 

O 

o 

o 

o 

O 

o 

o 

O 

O 

O 

ro 

to 

CD 

CM 

in 

go 

“ 

O 

ro 

CVi 

CM 

CM 

ro 

ro 

O 

8 


(S99JB9P)  010V  01  3AI1V33H  H10WIZV 


000*11  ooi'ot  0006  0008  OOOZ  0009  oo'os  000* 000  £ 


o 

ro 


o 

ID 

-L. 


o 

<7> 


o 

CJ 


o 

in 


o 

CD 


O 

CJ 


o 

o 

o 

O 

O 

o 

ro 

ID 

CJ 

CJ 

ro 

ro 

ro 

S31 


y- 


\ 


X 


X 


8S* 


O  ( 

<n  o 


X 


.§ 

8 


\ 


O 

O 

O 

O 

O 

o 

O 

O 

o 

p 

o 

O 

ro 

ID 

& 

CO 

in 

CD 

f- 

o 

ro 

CJ 

CJ 

CJ 

ro 

ro 

O 

8 


(saajBap)  QI0V  01  3AI1V13H  H10WIZV 


6000  7000  8000  9000  10,000  11,000 

DISTANCE  FROM  A0IO  (kilometers) 


o 

fO 


o 

tO 

JL 


o 
a > 


o 

CVJ 


O 

m 


o 

00 


o 

OJ 


o 

o 

o 

o 

O 

o 

fO 

to 

OJ 

OJ 

fO 

ro 

fO 

o 

o 

o- 


i 


o 

o_ 

o 

00 


ut  ^  ■ 
(N  O 


¥- 


■V 

t 

V, 


\ 


V 


X 


o 

-8 


o 

r 

o 

O 

o 

o 

o 

o 

6 

o 

270- 

o 

o 

to 

<o 

05 

0J 

in 

00 

0J 

CM 

o 

fO 

ro 

ro 

(/> 

k_ 

0) 

a) 

E 

o 


O 

o  § 
§  < 

O 

cr 

it. 

LU 

O 

sj 

8  52 

a 


8 


(sssjBap)  QI0V  01  SAIIVISH  HlflWIZV 


o 

s 


CN 

Q 

>i 

rd 

U 

u 

rd 

rQ 

3 

CO 


m 


a) 

Vh 

3 

tn 

rH 

P>4 


3  000  4000  5000  6000  7000  8000  9000  10,000  11,000 

—  1  1  ■ _ i— _ _ _ li _ t  i _ T  ■ 


O 

ro 


O 

CD 


O 

Oi 


O 

CJ 


o 

m 


o 

oo 


o 

CVJ 


o 

O 

O 

o 

O 

N 

o 

ro 

CD 

<\J 

OJ 

to 

to 

to 

o  o  o 


+ 


Y- 


\ 


V 


o  o  . 

00  ft  to  *-H 

HO  * 


X 


X 


o 

1-8 

O 


\ 


to 

k_ 

<D 

<D 

E 

o 


O 

o  § 
8  < 

o 

cc 


liJ 

o 

§< 
8  <2 
a 


o 

o 

O 

O 

O 

O 

O 

O 

o 

p 

O 

O 

ro 

CD 

0) 

CVJ 

1 n 

go 

2*. 

F» 

O 

to 

CM 

CVJ 

CVJ 

ro 

ro 

O 

8 


(saaj6ap)  QI0V  01  3AI1V33U  H10WIZV 


Q 

>i 

rd 

u 

u 

rd 

Xi 

D 

cn 


vD 


CD 

M 

3 

CP 

rH 

Cm 


3  000  4000 _ 5000  6000  7000 _ 8000 _ 9000  10,900  11,000 


o 

ro 


O 

to 


O 

0) 


O 

CVJ 


O 

to 


o 

00 


o 

CVJ 


o 

o 

o 

O 

O 

<r 

h* 

o 

ro 

to 

CJ 

CVJ 

ro 

ro 

ro 

y- 


\ 


V 


V 


o 

1-8 

O 


\ 


t/> 

w_ 

a> 

Q 

a) 

£ 

o 

>i 

rtf 

M 

H 

rtf 

o 

& 

<s» 

< 

P 

co 

s 

o 

cr 

r* 

U- 

r-H 

iu 

CD 

u 

U 

z 

P 

< 

tX> 

»- 

H 

CO 

Cl. 

Q 

o 

o 

O 

O 

o 

O 

o 

o 

o 

p 

o 

o 

ro 

to 

O) 

CVJ 

lO 

00 

F- 

o 

ro 

CVJ 

CVJ 

CVJ 

ro 

ro 

o 

8 


(S93J63P)  010V  01  3  Al  IV  3  3d  HlflWIZV 


3  0/30  4000  5000 _ 6000 _ 7000 _ 8000 _ 9000 _ IQ.pOO _ 11,000 


o 

ro 


O 

CO 


o 

O) 


o 

CVJ 


o 

to 


o 

00 


o 

CVJ 


o 

o 

O 

O 

O 

<fr 

o 

ro 

<0 

CVJ 

CVJ 

ro 

ro 

ro 

"T“ 

o 

ro 


~~r~ 

o 

CD 


y- 


\ 


v 


V 


V 


X 


\ 


O 

o 

O 

O 

6 

O 

■  T  " 

R 

CVJ 

o 

O 

o> 

CVJ 

in 

00 

CVJ 

cvj 

O 

ro 

ro 

ro 

(seej6ep)  QI0V  01  3AI1V13U  HlflWIZV 


</) 

k_ 

<D 

<D 

£ 

o 


O 

o  g 

O 

CC 

b_ 

LU 

O 


CO 

O 


w 

>1 

rtf 

rtf 

3 

CO 


00 
r— I 

0) 

u 

3 

0^ 


O 

s 


3  0,00  4000  5000  6000  7000  8000  9000  IOtQOO  11,000 


o 

O 

O 

o 

o 

o 

00 

o 

o 

o 

s 

o 

o 

O 

ro 

o 

CO 

o 

ro 

<0 

<n 

CM 

— 

CJ 

CM 

CM 

ro 

to 

to 

</> 

v- 

0) 

<D 

E 

o 


O 

<S> 

< 


O 

cr 


LU 

o 

< 

h- 

co 

o 


CM 

w 

>1 

03 

rQ 

P 

to 


cr» 


cu 

in 

P 

tn 

-H 

C*4 


3  000  4000  5000 _ 6000 _ 7000 _ 8000 _ 9000 _ 10,^00 _ 11 ,000 


o 

ro 


O 


O 

0> 


o 

CM 


o 


o 

00 


o 

cm 


o 

O 

o 

o 

O 

o 

ro 

CD 

CM 

CM 

ro 

ro 

ro 

"T“ 

O 

ro 


“T“ 

O 

CD 


y- 


\ 


X 


X 


X 


o 

1-8 

o' 


o 

■8 

0> 


\ 


6 

O 

O 

6 

o 

270- 

o 

O 

CM 

*o 

00 

CM 

M- 

CM 

o 

ro 

ro 

ro 

(/> 

k_ 

<D 

<1) 

E 

o 


O 

(3 


O 

cr 


LU 

a 

i? 

S  <2 

a 


nr) 

w 

>1 

nJ 

Vh 

u 

rQ 

3 

w 


o 

CN 

0) 

n 

3 

Cn 

H 

Cm 


o 

s 


(S99J69P)  010V  01  3AI1V138  HlfWIZV 


000*11  000*01  0006  0008  0001  0009  0005  000  OOoT 


O 

ro 


O 

<0 


O 

o> 


o 

OJ 


o 

m 


o 

oo 


o 

co 


o 

o 

o 

O 

O 

s 

o 

ro 

<0 

OJ 

00 

ro 

ro 

ro 

“T“ 

O 

ro 


4 


y- 


x 


X 


v 


X 


\ 


o 

to 


o 

0> 


o 

00 


o 

in 


o 

oo 


o 

00 


O 

P 

O 

O 

*3“ 

F- 

O 

ro 

0J 

co 

ro 

ro 

O 

g 


(S99J69P)  010V  01  3AI1V13U  Hinwizv 


5000  6000  7000  8000  9000  10,000  lljOOO 

DISTANCE  FROM  A0IO  (kilometers) 


3  000  4000  5000  6000  7000  8000  9000  10,000  11,000 

■  1.  I  i  — —  ■  ...» i _  I  i  _ 7  1 - 


O 

ro 


O  o 

<0  cn 

1  1 


o 

CM 


O 

m 


o 

CO 


o 

CM 


o 

O 

O 

O 

O 

M- 

N 

o 

ro 

w 

CM 

CM 

ro 

ro 

ro 

n  o 


~~r~ 

o 

ro 


O 

<£> 


O 

<T> 


O 

CM 


y- 


\ 


V 


V 


-?4- 


o 

ID 


O 

00 


o 

CM 


X 


o 

■8 


\ 


O 

■8 

0> 


w 

k_ 

CD 

rH 

+- 

0> 

Eu 

£ 

_o 

>1 

rtf 

-ie 

M 

W' 

M 

O 

rtf 

G» 

3 

< 

cn 

O 

cr 

u_ 

CM 

CM 

UJ 

u 

Q) 

U 

< 

3 

c n 

rH 

a 

o 

O 

O 

O 

M- 

r- 

O 

ro 

CM 

CM 

ro 

ro 

O 

s 


(S99J69P)  010V  01  3  AI1V13H  HlflWIZV 


3  000  4000  5000  6000  7000  8000  9000  10,000  11,000 


o 

ro 


O 

to 


O 

CD 

_L_ 


O 

CVJ 


O 


o 

oo 


o 

CVJ 


o 

o 

O 

O 

O 

<*■ 

N 

o 

ro 

to 

CVJ 

CM 

ro 

ro 

ro 

2  s : 


-A 


y- 


\ 


v 


V, 


s 


58' 


X 


o 

1-8 

o' 


\ 


to 

a> 

CM 

to 

Cm 

E 

o 

X 

rtJ 

M 

w 

M 

O 

rd 

<S> 

rQ 

< 

3 

cn 

O 

cr 

u_ 

ry 

r\ J 

UJ 

a 

(1) 

z 

< 

3 

h- 

t7» 

co 

Q 

Cm 

o 

h8 

to 


o 

o 

O 

o 

o 

O 

o 

O 

O 

R 

O 

O 

ro 

to 

CD 

CVJ 

«n 

CD 

CVJ 

<r 

O 

ro 

cvj 

cvj 

ro 

ro 

o 

s 


(S33J63P)  010V  01  3  AllV13d  Hinwizv 


000‘n  OOOOI  0006  0008  000 1  0009  OOOS  000*  000  £ 


o 

O 

o 

o 

o 

o 

CO 

o 

o 

o 

s 

o 

o 

o 

ro 

o 

<0 

o 

ro 

<0 

0) 

OJ 

— 

OJ 

CO 

OJ 

ro 

ro 

ro 

</> 

w. 

CL> 

•*— 

<D 

E 

o 


o 

< 

o 

IT 

Ll 

UJ 

o 


CO 

Q 


(saaj6ap)  QI0V  01  3AI1V13H  HlflWIZV 


on 

>1 

m 

CO 


OJ 

CD 

3 

tr 

rH 

C^| 


3  000  4000  5000  6000  7000  8000  9000  10,000  11,000 


o 

lO 


o 

CD 


O 

O) 


O 

CM 


o 

m 


o 

00 


o 

CJ 


o 

o 

o 

o 

O 

X 

o 

ro 

CD 

CM 

CM 

rO 

ro 

ro 

V- 


X 


X 


X. 


X 


\ 


O 

O 

O 

O 

O 

o 

o 

O 

o 

D 

o 

O 

ro 

CD 

<T> 

C\J 

m 

go 

— 

rf 

X 

O 

ro 

CM 

CM 

CM 

ro 

ro 

(S99J69P)  010V  01  3 AI1V13H  HinWIZV 


o 

-8 


Q 

W 

k_ 

g 

<D 

00 

CD 

6 

o 

>1 

rtf 

u 

u 

rtf 

O 

o 

8 

5) 

< 

3 

v> 

X 

o 

cr 

LO 

Ll. 

OO 

UJ 

0) 

o 

8 

< 

H 

H 

O 

CD 

C/) 

tu 

O 

-8 

K> 


O 

8 


APPENDIX  A 


SEISMIC  DATA  LABORATORY 
ALEXANDRIA,  VIRGINIA 


DIGITAL  COMPUTING  SECTION 


A .  IDENTIFICATION 
Title:  ANOMALY 

COOP  Identification:  Z108 
Category.  General 
Programmer:  Nicholas  Fletcher 

Date:  5  July  1967 

B „  PURPOSE 

To  compute  travel-time  and  amplitude  anomalies 
C .  USAGE 

1.  Operational  Procedure:  This  is  a  FORTRAN  63  program.  All 
input  is  from  punched  cards  and  the  output  is  from  the  printer, 
and  an  off-line  card  image  tape 

2.  Parameters :  Input  parameters  consist  of  the  latitude  and  longi¬ 
tude  of  the  location  of  the  subarray  centers. 

3  Space  Reguired:  Not  Applicable 

4  Temporary  Storage  Reguirement s :  None 

5.  Alarms .  None 

6.  Error  Codes:  None 

7.  Error  Stops:  None 

8  Input  and  Output  Tape  Mounting.  There  is  no  input  tape.  The 
output  tape  is  on  logical  Unit  16. 

9 .  Input  and  Output  Formats: 


-  a.: 


The  input  card  formats  are  as  follows. 

Subarray  locations  deck,  must  be  the  first  deck 
of  cards  read  by  the  program.  The  first  card  of 
this  deck  must  contain  the  number  of  subarray 
location  cards.  This  number  is  punched  in  Columns 
1  through  5  right  justified.  The  first  card  is 
followed  by  the  subarray  location  cards,  with  one 
location  per  card.  The  program  is  limited  to  a 
maximum  of  51  subarrays,  1  (One)  punched  in  Column 
10  prints  out  J.  B„  Table  and  station  locations. 


Subarray  location  card  format  is  a_s  .follows; 


FORMAT  (A8,  2  (F5,0  F3  0  F5.1,  Al) 

COLS  „ 


1-8  Subarray  name,  left  gustified 

9-13  Subarray  latitude  degrees 

14-16  Subarray  latitude  minutes 

17-21  Subarray  latitude  seconds  (with  decimal  point 

punched ) 

22  N  or  S  for  north  or  south  latitude 

23-27  Subarray  longitude  degrees 

28-30  Subarray  longitude  minutes 

31-35  Subarray  longitude  seconds  (with  decimal  point 

punched ) 


36  E  or  W  for  east  or  west  longitude 

The  location  cards  are  followed  by  the  first 

event  card  - - 

Event  card  format  is  as  follows . 

Columns  1  through  14  contain  the  origin  time  of  the 
event „ 

C_OLS, 

1~2  Month 

3-4  Day 


-  A  -  2 


COLS « 


5-6  Year 

7-8  Hour 

9-10  Minute 

11-14  Second  to  the  nearest  tenth  of  a  second 

15-18  Latitude  of  the  event  in  degrees  to  the  nearest  tenth 

19  N  or  S  for  North  or  South 

20-24  Longitude  of  the  event  in  degrees  to  the  nearest  tenth 
25  E  or  W  for  East  or  West  longitude 

26-28  Depth  of  event  to  the  nearest  kilometer 
29-31  Magnitude  to  the  nearest  tenth 

32-36  This  is  an  output  field.  The  program  will  compute 
the  distance  (to  the  nearest  kilometer)  between 
the  event  and  the  standard  subarray  (AO  for  LASA) 

37-39  This  is  an  output  field.  The  program  will  compute 
the  station  to  epicenter  azimuth  (to  the  nearest 
degree  relative  to  AO) 

40-78  Not  used  by  the  program,  but  reproduced  on  output. 

79-80  Must  be  blank,  to  identify  this  card  as  an  event 
card  and  not  a  phase  card . 

Each  event  card  should  be  followed  by  phase  cards. 

A  phase  card  contains  the  observed  phase  data  for 
a  given  subarray.  Each  phase  card  must  have  an 
index  number  punched  in  Columns  79-80.  This  index 
is  established  by  the  order  of  the  subarray  location 
cards  (see  Sec-C9) .  There  does  not  have  to  be  a 
phase  card  for  every  subarray  location  card  nor  a 
specific  numerical  order. 

Phase  Card  Format : 

FORMAT  (Columns  1  through  9  contain  the  phase  arrival 
time  at  this  subarray) 


-  A  -3 


Hour 


COLS  . 
1-2 
3-4 
5-9 
10 

11-14 

15-16 

17-20 

21 

30 

31-35 

36-40 

41-42 

43 

44-45 

46-50 

51-55 

56-60 

61-78 

79-80 


Minute 

Second  to  the  nearest  hundredth 
Not  used 

Amplitude  (zero  to  peak)  to  the  nearest  millimicron 
Not  used 

Amplitude  (peak  to  peak)  to  the  nearest  millimicron 
Not  used 
Signal  Quality 

Time  anomaly  to  the  nearest  hundredth  of  a  second  — 
this  is  an  output  item  computed  by  the  program. 

Amplitude  anomaly  zero  to  peak  relative  to  geometric 
mean  computed  by  the  program  (output  field) 

Subarray  name  code  (output  item) 

Not  used 

Reference  subarray  (for  each  event)  (only  in  last 
phase  card) 

Peak  to  peak  amplitude  anomaly  relative  to  geometric 
mean  (output  field) 

Zero  to  peak  amplitude  anomaly  normalized  to  average 
amplitude  (0-P)  (output  field) 

Peak  to  peak  amplitude  anomaly  normalized  to  average 
amplitude  (P-P)  (output  field) 

Reserved  but  not  used  by  program 

Index  number  of  subarray  recording  this  phase  must 
correspond  to  subarray  location  card. 

Additional  events  with  their  associated  phase  cards 
follow  in  like  fashion  for  as  many  events  not  repeated 
in  the  deck  setup. 


-  A  -4 


The  output  tape  is  a  card  image  tapef  with  the  same 
information  and  format  as  the  input  tape  with  the 
addition  of  the  computed  values  for  distance  and 
azimuth  for  event  cards,  and  travel-time  and  ampli¬ 
tude  anomalies  for  the  phase  cards, 

10.  Selective  Jump  and  Stop  Settings;  None 

11.  Timing :  Approximately  3.5  seconds  per  event 

12 .  Accuracy.  No  better  than  the  input  data 

13.  Cautions  to  User:  Every  phase  card  must  have  an  index 
number.  If  Columns  79-80  are  blank,  the  program  assumes 
it  is  an  event  card. 

14.  Equipment  Configuration:  Standard  F-63  COOP  monitor 
system  with  the  card  reader  on  Unit  7  and  the  printer  on 
Unit  6  and  the  output  tape  on  Unit  16. 

15.  References :  Letter  from  MIT  dated  April  8,  1966 

D .  METHOD 

Def  mit  ions 

Observed  Travel-Time  =  Observed  Phase  Arrival  Time  Minus  Event 

Origin  Time 

Predicted  Travel-Time  =  Travel-Time  for  the  J-B  Table.  Found 
by  table  look-up  on  distance  and  depth  with  four-point 
interpolation  for  distance  and  linear  interpolation  for 
depth . 

0  =  Observed  travel-time  at  subarray  i 

H.  =  Predicted  travel-time  to  subarray  l 

l 

O  =  Observed  travel-time  at  standard  subarray 
r 

H  =  Predicted  travel-time  to  standard  subarray 
r 

Transit  Time  Residual  (Travel-Time  Anomaly)  =0  -0  +H  -H, 

i  r  r  i 

(Relative)  Amplitude  Anomaly  =  Amplitude  (Zero  to  Peak)  at  sub¬ 
array  l  divided  by  geometric  mean  amplitude  (zero  to  peak) 
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(Average)  Amplitude  Anomaly  =  Amplitude  (Zero  to  Peak)  at 

subarray  1  divided  by  arithmetic  average  amplitude  (zero 
to  peak)  for  all  functioning  subarrays. 

Relative  and  average  amplitude  anomalies  are  also  computed  in 
like  manner  for  the  peak  to  peak  measurements 
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APPENDIX  B 


SEISMIC  DATA  LABORATORY 
ALEXANDRIA ,  VIRGINIA 

DIGITAL  COMPUTING  SECTION 


A .  IDENTIF  ICATION 

Title;  Display  and  Analysis  of  Anomalies 
COOP  Identification:  Z95  DISPLAY 
Category;  General 
Programmer:  J.  W.  Monroe 

Date :  19  July  1966 

B .  PURPOSE 

To  take  the  output  from  the  anomaly  program  and  compute 
the  average,  standard,  deviation  and  number  of  occurences  of 
a  selected  anomaly  for  selected  subarrays  of  sensors,  distances, 
and  azimuths „ 

C „  USAGE 

1.  Operational  Procedure  This  FORTRAN-63  program  is  with  a 

binary  or  symbolic  deck  and  a  set  of  control  cards  for 
each  use 

2.  Parameters.  None 

3.  Space  Reguired:  21  222  buffers 

4.  Temporary  Storage  Reguirement:  None 

5.  Prmt-Out s :  The  program  first  prints  out  the  anomaly  used, 
followed  by  the  distance  and  azimuth  criteria  and  the  sub¬ 
arrays  or  sensors  used  If  any  seismograms  are  to  be  deleted 
from  processing,  they  are  listed  next.  Following  this  will 
be  a  page  for  each  subarray  or  sensor  with  the  average 
standard  deviation  and  number  of  occurrences  of  the  selected 
anomaly  for  each  distance  and  azimuth  chosen.  If  plotting 

is  requested,  the  last  page  will  contain  scale  factors  and 
other  information  pertaining  to  the  plots „ 
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Error  Returns:  None 


6 . 

7 .  Error  Stops:  None 

8 .  Input  and  Output  Tape  Mountings. 

Logical  Unit  #4  -  Input  BCD  tape  from  ANOMALY 
Logical  Unit  #3  -  Output  plot  tape  if  desired 


9. 


Input  and  Output  Formats: 


INPUT: 

(only 

one  case 

per  run) 

Card  # 

Cols 

Format 

Description 

1 

1 

11 

Anomaly  to  be  used  (see  last 

page) 

5 

11 

Plot  switch  1-plot,  blank-no 

plot 

2 

1-10 

110 

Number  of  Dels  (Max  19) 

11-20 

110 

Number  of  Azimuths  (Max  13) 

21-30 

110 

Number  of  Subarrays  or  Sensors 

(Max  21) 

next  N 

1-10 

F10.2 

Distance  (km) 

11-20 

F10.2 

+_  Distance  (km) 

Where  N 

=  No. 

punched  in  Cols  1-10  on  card  #2  (Max  19) 

next  M 

1-10 

F10.2 

Azimuth  (Deg) 

11-20 

F10.2 

+  Azimuth  (Deg) 

Where  M 

=  NO. 

punched 

in  Cols  11-20  on  card  #2  (Max  13) 

next  L 

1-10 

FlOo  0 

Subarray  or  Sensor  Number 

Where  L 

=  NO. 

punched 

in  Cols  21-30  on  card  #2  (Max  21) 

next  D 

1-4 

14 

Seismogram  Nos.  to  be  deleted 

(Max  100) 


OUTPUT :  The  printed  output  will  be  that  described  in 

PRINT-OUTS  The  plot  tape  on  logical  unit  #3  will  con¬ 
tain  a  plot  for  each  subarray  or  sensor  of  the  average 
for  that  subarray  minus  the  average  of  averages  for  all 
subarrays.  The  plotted  value  will  be  computed  for  each 
distance  at  each  azimuth  (in  that  order)  until  N  points 
have  been  plotted  where  N  equals  the  number  of  dels  times 
the  number  of  azimuths.  The  plots  will  always  represent 
a  new  del  every  tenth  of  an  inch. 
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10.  Selective  Jump  and  Stop  Settings:  None 

11.  Timing :  2  mili-sec  per  case  where  each  case  has  a  separate 

distance,  azimuth  subarray,  and  anomaly. 

12 .  Accuracy:  The  same  as  input  data 

13.  Cautions  to  Users.  Both  logical  units  3  and  4  must  be 

mentioned  in  the  COOP  card.  If  no  plot  is  desired, 
use  the  BY3  convention.  The  plot  routine  is  not 
internally  controlled  so  the  user  must  request  the 
same  number  of  files  to  be  plotted  as  there  are  subarrays 
or  sensors  with  a  maximum  of  ten  across  the  plot .  Observe 
that  there  are  maximum  limits  on  the  number  of  distance 
and  azimuth  intervals  and  number  of  subarrays  and  seis¬ 
mograms  to  be  deleted . 

14.  Equipment  Configuration:  Standard  COOP  for  FORTRAN-63 

15.  References: 

Program  write-up  for  ANOMALY  Fletcher  -  16  June  1966 


METHOD 


Average  =X.  .  .  =  -  )  X 

l ,  3  ,k  N  .  L  l ,  3  ,k 


Standard  Deviation  =  a. 

i ,  3  ,k 


r  V*  2  ^  f  \2  ,  -I 

N  -I 

1*3 


Number  of  Occurrences  =  N 

i  *  3  *k 


where,  l 

3 

k 


azimuth  +_  deg's, 
distance  +  km's 
subarray  or  sensor 


N  =  number  events  pertinent  to  the  specified  coordinates 
and  specified  subarray 


X. 

i*3 


*k 


=  anomaly  used 
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MEMORANDUM 


TO: 

FROM. 

DATE: 

SUBJECT: 

REFERENCE: 


Users  of  Program  DISPLAY 
J.  Wo  Monroe 
2  September  1966 

Additional  Features  of  Program  DISPLAY 
Z95  DISPLAY  Program  Writeup  19  July  1966 


The  user  may  now  run  DISPLAY  using  more  than  one  input  tape  de¬ 
scribed  in  the  program  writeup  referenced  above  The  first  input  tape 
must  be  on  Logical  Unit  4  but  the  second,  third  etc. ,  tapes  may  be 
assigned  to  any  logical  unit  from  1  thru  49  inclusive,  except  Logical 
Unit  3  (the  plot  tape) .  When  using  more  than  one  input  tape.  Logical 
Units  3,4,  and  all  other  logical  units  must  be  mentioned  in  the  COOP 
card.  The  control  card  setup  for  the  case  of  more  than  one  input  tape 
is  exactly  the  same  as  that  shown  in  the  SAMPLE  INPUT  DATA  FOR  DISPLAY 
of  the  writeup,  except  that  following  the  last  card  of  seismograms  to 
be  deleted  or  following  the  last  subarray  number  if  there  are  no  de¬ 
letion  cards,  a  blank  card  must  be  inserted  followed  by  the  logical 
tape  assignment  for  the  second  input  tape.  If  more  than  two  input 
tapes  are  used,  additional  cards  are  used  to  specify  the  logical  tape 
unit  for  input  tapes  3,4,  etc.,  but  the  blank  card  is  not  used  again. 
Format  for  these  tape  assignments  is  15.  The  following  is  a  sample 
card  setup  starting  with  the  deletion  cards. 


etc . 


fourth  input  tape 


logical  tape  unit  for  the 
third  input  tape 

logical  tape  unit  for  the 
second  input  tape 


signifies  the  end  of  the 
^seismograms  to  be  deleted 
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As  soon  as  an  end  of  file  is  encountered  on  the  input  taps  being  pro¬ 
cessed,  a  tape  assignment  card  is  read  in  and  the  program  pauses  with 
the  next  input  logical  unit  number  displayed  in  the  "A  register.  When 
the  Start  Key  is  hxt ,  the  program  begins  processing  the  new  input  tape. 
Two  input  tape  units  might  be  used,  such  as  4  and  5  where  the  first 
tape  is  on  Unit  4  and  the  second  tape  is  mounted  on  5  while  4  is  being 
processed.  After  Unit  5  begins  processing,  another  could  be  mounted 
on  Unit  4. 


DEBUGGING  AIDS  FOR  DISPLAY 


Jump  Key  2 


Jump  Key  1 
Jump  Key  3 


On  displays  the  event  number  being  processed  in  octal 
in  the  "A1  register.  Program  halts  until  Start  Key 
is  hit.  For  exact  stops,  leave  JK-2  cn  and  step 
through  events  with  Start  Key.  Program  halts  on  header 
card  of  event  being  displayed  in  A'  register. 

Up  causes  printout  of  accumulation  that  far  into  the 
program.  Turn  JK-1  off  as  soon  as  printout  begins 
and  JK-3  off  as  soon  as  printout  has  finished.  Use 
this  feature  sparingly  since  some  loss  of  accuracy 
occurs  when  sigma  is  set  equal  to  zero  rather  than  a 
negative  number.  All  other  times  the  accumulation 
is  reset  to  its  original  value  without  loss  of  accu¬ 
racy  . 


Be  sure  both  COOP  card  and  EXECUTE  card  have  sufficient  time  limits  for 
the  case  being  run.  For  a  case  of  324  events  with  21  stations  17 
distances,  and  12  azimuths  allow  at  least  60  minutes 


CHANGE  IN  FORMAT 


The 
grams  to 


format  of  the  control  cards  containing 
be  deleted  has  been  changed  from  14  to 


the 

15 


numbers  of 


seismo- 
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